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Abstract 

Hamiltonian mechanics can be used to constrain temperature simultaneously with energy. We 
illustrate the interesting situations that develop when two different temperatures are imposed 
within a composite Hamiltonian system. The model systems we treat are (f> chains, with quartic 
tethers and quadratic nearest-neighbor Hooke's-law interactions. This model is known to satisfy 
Fourier's law. Our prototypical problem sandwiches a Newtonian subsystem between hot and cold 
Hamiltonian reservoir regions. We have characterized four different Hamiltonian reservoir types. 
There is no tendency for any of these two-temperature Hamiltonian simulations to transfer heat 
from the hot to the cold degrees of freedom. Evidently steady heat flow simulations require energy 
sources and sinks, and are therefore incompatible with Hamiltonian mechanics. 
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INTRODUCTION 



Molecular dynamics has facilitated the modeling of the macroscopic equation of state and 
constitutive relations in terms of atomistic Hamiltonians. Three research groups developed 

nn 

the method[lH3( and its applications: Fermi's one- and two-dimensional simulations at Los 
Alamos, Vineyard's radiation-damage simulations at Brookhaven, and Alder's if -Theorem 
and phase-change studies at Livermore. All three groups treated the approach to equilibrium 
from nonequilibrium initial conditions. Such problems are straightforward applications of 
Newtonian, or Lagrangian, or Hamiltonian mechanics. Today both the scale and the com- 
plexity of the modeling have expanded to make molecular dynamics a comprehensive tool 
for learning and understanding. 

Some conceptual problems remain|4|. Irreversibility, and its connection to Lyapunov 
instability - the exponential growth of small perturbations - are undergoing wide- 
ranging investigations. Nonequilibrium boundary conditions have been under intensive 



development 
definition o 

temperature 



From a theoretical standpoint there is no agreement on the proper 
nonequilibrium states. How should they be described? Nonequilibrium 
^ can appear to be a difficulty The many equivalent equilibrium definitions, 
thermodynamic, kinetic, and configurational, all differ from one another away from equilib- 
rium, with the latter two temperatures becoming tensors rather than scalar quant it iesflOf] . In 
strong shockwaves the longitudinal and transverse kinetic temperatures can differ temporar- 
ily by an order of magnitude HI] , but then equilibrate in a few collision times. Sometimes 
the local configurational temperature is negative 12] . 

Stationary nonequilibrium flows require nonequilibrium boundary conditions to impose 
local velocities and temperatures and to extract the irreversible heat generated by equili- 
bration pr o cesses {5, 0]. The Nose-Hoover version of Nose's constant-temperature dynamics 
provides a robust time-reversible approach to temperature control. That approach uses inte- 



gral feedback. Isokinetic velocity rescaling 13j can likewise be described with time-reversible 



deterministic motion equations using differential feedback. Couette shear flow, Fourier heat 
conduction, and shockwave propagation are prototypical examples of nonequilibrium prob- 
lems which can all be driven and maintained by special boundary regions with stationary 
velocities and temperatures 4[] . 

Because there are many definitions of temperature it is natural to explore their relative 



2 



usefulness in driving flows away from equilibrium. Both kinetic and configurational tem- 
peratures can be controlled by constraining Hamiltonian systems. Here we c omp are five 
such approaches to the simplest heat-flow problem, conduction in a 4 chain 14j-|l7j . Such a 
chain combines harmonic Hooke's-Law interactions with quartic tethering potentials. The 
4 chain is an improved relative of the seminal anharmonic chain models studied by Fermi's 
group at Los Alamos. His models lacked the realism of the tethered-particle simulations 
studied in much greater detail fifty years later, by Aoki and Kusnezov. Just as Fermi's 
results (lack of equilibration) were a surprise to Fermi, the present results (once again, lack 
of equilibration!) surprised us. What we find is apparently a common and serious drawback 
of Hamiltonian thermostats, a failure to promote heat flow 171. The failure of recent workers 

n 

to recognize this defect [18] helped motivate the present work. 

The body of this paper has three parts. In Section II we describe the thermostats to 
be considered. In Section III we compare the results of sandwich simulations of [ cold + 
Newtonian + hot ] thermostated systems. Section IV is a summary of our findings and the 
conclusions which we draw from them. 



SEVERAL HAMILTONIAN THERMOSTATS 



Nose|19l-l21| took a bold step forward in 1983, finding a Hamiltonian mechanics which 
could model the canonical-ensemble isothermal distribution. His Hamiltonian, for # degrees 
of freedom , among which is a new "time-scaling" variable s (we use s 2 here to ensure that 
the logarithm is meaningful) along with its conjugate momentum p s , is : 

•Hnosc = Y> 2 /2ms 2 ) + Hq) + (p 2 s /2M) + (#fcT/2) ln( S 2 ) . 

The adjustable parameter M can be used to vary the timescale of the ( s,p s ) thermostat. 
Nose showed that the dynamics for this ( q,p,s,p s ) Hamiltonian was consistent with the 
stationary Gibbsian canonical distribution for the coordinates and the "scaled momenta" 
{ q,(p/s) } : 

^Nosc — > f(q,p,s) ex exp[ -(K(p/s)/kT) - (Hq)/kT) ] . 



A simpler route [22] to a similar result is to begin with the Nose- Hoover motion equations, 
a modification of Nose's work in which s is completely absent and p s is replaced by a time- 
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reversible friction coefficient, (, along with its characteristic relaxation time r : 
{ q = (p/ m ) ; p = F(q) - (p } ; ( = (p 2 /mkT) - 1 ]/r 2 . 

Next, it is easy to see that a Gaussian distribution for (, along with the canonical distribution 
for the (q,p) variables satisfies the phase-space continuity equation (a generalized Liouville 
Theorem), giving the result : 

(df/dt) = — ► f{q,p,() oc exp[ -(K{p)/kT) - (<S>(q)/kT) - (l/2)(Cr) 2 ] . 

£ controls the flow of energy to and from the thermostated system with the arbitrary relax- 
ation time t. Bauer, Berry, Brahka, Bulgac, Hamilton, Jellinek, Klein, Kusnezov, Martyna, 



Tuckerman, Wink 
tion equations 



2J29 



er, and Wojciechowski all suggested various modifications of these mo- 



29|. Surprisingly, Bulgac and Kusnezov were able to demonstrate the 



feasibility of modeling Brownian Motion with strictly time-reversible deterministic motion 
equations by introducing three thermostat variables like ( , rather than just one or two 30]. 



Lagrangian and Hamiltonian mechanics suggest a wide variety of approaches to ther- 



mostating. Dettmann and Morriss 



31 



32] discovered Hamiltonian bases for both the Nose- 



Hoover and Gaussian isokinetic motion equations. Bond, Laird, and Leimkuhler rediscovered 
this approach a year later {33)]. [ The isokinetic equations are the instantaneous r — > limit 
of the Nose-Hoover motion equations given above. ] Landau and Lifshitz' "configurational 
temperature" 34j was rediscovered and generalized by Braga, Rugh, and Travis [7|, |8|, |35| , 



kT c = (F 2 )/( V 2 U ) . 

Constant configurational temperature can be imposed as a straightforward ( and tedious ) 
holonomic constraint (3, S H] ■ 

Recently Campisi, Hanggi, Talkner, and Zhan suggested the use of a logarithmic 
thermostat 361] . very much like Nose's, but without an explicit coupling to the remaining 
degrees of freedom : 

•Hczth = (P 2 s /2m) + (kT/2) \n(s 2 + 5 2 ) . 

Here we include a S 2 in the definition so as to avoid divergence when s changes sign. The 
equations of motion of the logarithmic thermostat, absent the coupling forces linking it to 
the system it influences, are : 



s = ip s /m) ; p s = -skT/(s 2 + 5 2 ) 



(sp s ) = ( (d/dt)(s Ps ) )-(sp a ) = -(sp a )~ -kT . 

Notice that the time-averaged derivative of a bounded quantity, here (sp s ), vanishes. This 
averaging operation shows that when 5 can be ignored the time-averaged kinetic energy 
of the "thermostat" is ( (p 2 s /2m) ) = (kT/2) . The relatively poor results obtained with 

-spanol, Hoover, and Melendez and documented in 
. It is worth noting here that the instantaneous 



this thermostat were pointed out b 
a series of ar^iv contributions [isl, ElZ 
configurational temperature of the pure logarithmic-thermostat potential is negative, — kT ! 

<j){s) = kTln(s) -> { F(s) = (-kT/s) ; V 2 = (-kT/s 2 ) } ->■ kT c = -kT ! 

This bipolar character, with T c = — ( T& ), is definitely uninviting. 

An isokinetic form of mechanics, based on a nonholonomic ( velocity-dependent ) con- 
straint, can be imposed by using the Hoover-Leete Lagrangian and Hamiltonian 
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£hl( Q,<j) K(q) = ^(mg 2 /2) constant 



^hl( q,p) = 2yjK(p)K(q) - K(q) + $(g) . 

This Lagrangian approach closely resembles the Gauss'-Principle isokinetic approach, though 
the resulting trajectories are quite different. In the Hoover-Leete approach the momenta 
{ p } evolve in the usual way, but the velocities are continuously rescaled, by a Lagrange 
multiplier A, so that there are two different versions of the kinetic energy : 

{ 4 = (p/m)/(l + A) ; p = F } ; K(p) = ]>> 2 /2m) ^ £(™f/2) = K(q) . 

In the remainder of this work we apply the original Nose-Hoover thermostat, Nose's 
thermostat, the logarithmic thermostat, and the Hoover-Leete thermostat, to a simple model 
system known to follow Fourier's law, the 4 model, so called because each particle in a 
harmonic chain is, in addition, tethered to its lattice site by a quartic potentia. 
Kusnezov carefully characterized the 4 model's dependence on temperature fl4 
the masses and both force constants and Boltzmann's constant all set equal to unity the 4 
model has a heat conductivity k ~ 2.8/T 4 / 3 in one dimension (l5l]. 

COMPARING FIVE CONDUCTIVITY APPROACHES IN <h 4 CHAINS 



Aoki and 
Q. With 



The 4 model can be implemented in any number of dimensions and with any lattice 
structure. Likewise the number of particles and the size of the temperature gradient can be 
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large or small. Fortunately Aoki and Kusnezov have carried out a series of comprehensive 
investigations establishing that the model is a useful representation of Fourier conductivity 
over a wide range of conditions and dimensionality. The temperature profiles they found 
have a typical "jump" at system boundaries but are otherwise quite unremarkable [14. [15]. 
For simplicity we confine all of our investigations reported here to a 60 particle chain, with 
twenty Newtonian particles sandwiched between a cold 20-particle reservoir, at temperature 
0.5 and a hot 20-particle reservoir at a temperature 1.5. We choose an overall energy such 
that the Newtonian particles start out with an average temperature of order 1.0. Except 
for the Nose-Hoover case, all of the simulations carried out here have constant total energy, 
due to their Hamiltonian character. 



Nose-Hoover Thermostat 



The straightforward approach to this problem uses the ( nonHamiltonian ?) Nose- Hoover 



equations of motion 2l|, |22| in both 20-particle reservoir regions : 



{ q = (p/m) ; p = F{q) - (p } ; C = [ ( ( P 2 /mkT) ) - 1 ]/r 2 . 

In our numerical work the particles have unit mass. Also, both Boltzmann's constant k and 
the thermostat relaxation times ( r co id, Th t ) have been chosen equal to unity. We indicate 



the questionable nonHamiltonian character here, in order to emphasize that Dettmann|31| 
(and later, Dettmann with Morriss[32|) found and described a special Hamiltonian which, 
when set equal to zero, generates these same equations of motion, but only for a single 
temperature — not simultaneously for two or more of them. 

At any time t the average values of p 2 for the two reservoir regions, indicated here by 
( p 2 ) , are one-twentieth the sum of the 20 instantaneous cold or hot reservoir-particle 
contributions. The forces in this problem are of two kinds: nearest-neighbor Hooke's-Law 
forces and cubic forces from the quartic tethering potential : 

$ = £(l/2)(x i+1 - 1.0 - Xi ) 2 + £(l/4)(x, - x m f , 

where the fixed sites { x Q } have a regular lattice spacing of unity. Likewise the particle 
masses, force constants, relaxation times, and Boltzmann's constant have all been set equal 
to unity, for simplicity. 
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FIG. 1: Kinetic and configurational temperature profiles for a 60-particle (f) A chain according 

to the Nose-Hoover equations of motion. The particle temperatures are averages over a billion- 
timestep simulation with dt = 0.01. By "billion" we use the word to mean one thousand million, 
10 9 , throughout this work. Here, as in a^of our simulations, we use the classic fourth-order Runge- 
Kutta integrator. The time-averaged heat flux is -0.0618, with the flow from right to left. This 
corresponds to a heat conductivity k ~ 2. Simulations with longer chains are in agreement with 
Aoki and Kusnezov's work cited in References 14 and 15. 

Figure 1 shows the kinetic and configurational temperature profiles that result from a 
billion-timestep fourth-order Runge-Kutta solution of the 20 "cold" plus 20 "Newton" plus 
20 "hot" pairs of equations for { q,p } plus the additional two differential equations for the 
cold and hot friction coefficients { Ccoid, Chot } • The resulting kinetic temperature profile is 
unremarkable, save for the slight jumps at the system boundaries and thermostat interfaces. 
The difference between the two temperature definitions, kinetic and configurational, is a 
measure of the need for clarity in defining nonequilibrium system properties. Nose-Hoover 
integral feedback guarantees precise temperature control. The kinetic temperatures are 
automatically equal to the target temperatures ( 0.5 in the cold region - 1.5 in the hot 
). This follows because the long time averages of the differential equations for the friction 
coefficients necessarily reproduce exactly the specified kinetic temperatures : 
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FIG. 2: Kinetic and configurational temperature profiles using Nose's original Hamiltonian without 
time scaling to constrain the temperatures of 20 "Cold" and 20 "Hot" particles. In each of the 
three regions the configurational temperatures (open circles) are in good agreement with the kinetic 
temperatures. This simulation fails to generate a heat flux from hot to cold. The time averages 
come from a billion-timestep simulation with alt = 0.0005. The value of the Hamiltonian is 15.000. 

Carlos Braga and Karl Travis developed a similar automatic feedback approach for con- 
figurational temperaturej^. The heat flux for this problem corresponds to a conductivity 
significantly smaller than the large system limit, k ~ 2.8, established by Aoki and Kusnezov. 
We have confirmed that longer chains agree very well with Aoki and Kusnezov's work. Let 
us use this 60-particle solution as a "standard case" for the 4 model and investigate the 
same problem using four types of Hamiltonian heat reservoirs. 



Original Nose Thermostat (with s included in the Motion Equations) 



In Nose's revolutionary 1984 publications ljj, |20) he introduced his original Hamiltonian, 
which contains an additional conjugate pair of Hamiltonian variables (s,p s ) : 



H No se = J2(P 2 / 2ms2 ) + + (#^7 2 ) ln(s 2 ) + (p*/2M) 



Again choosing the path of simplicity, we set the degrees-of-freedom parameter and the 
thermostat's effective mass M both equal to 20. We set the particle mass m and Boltzmann's 
constant both equal to unity. The results which we find for Nose's original Hamiltonian 
thermostats are shown in Figure 2. The relative stiffness of the motion equations for the 
particles , 

# 

{ q = (p/ms 2 ) ; p = F } ; s = (p s /M) ; p s = (p 2 /ms 2 ) -kT]/s , 

is a consequence of Nose's time-scaling variable s , which appears in the denominator. The 
stiffness requires a much smaller timestep dt. Problems with only a few degrees of freedom 
tend to be singular unless (s 2 ) is replaced by (s 2 + 5 2 ), a precaution not needed here. We 
used dt = 0.0005 in order to conserve energy to six-digit accuracy over the course of a 
billion-timestep run. 

Despite the extra work due to the smaller timestep dt, twenty times more for the same run 
duration, relative to the Nose-Hoover equations, the temperature "profile" is disappointing. 
It shows a lack of effective interaction between the thermostated regions and the Newtonian 
Partlcl e, These results are tyP M- Higher or .ower energ.es, or longer or 8 ho rt er s y ,e m8 , 
provide similar results, establishing that Nose's Hamiltonian is ineffective for heat flow 
problems. The two "time-scaling" variables ( s co id, Shot ), with initial values of unity, had 
averaged values of 16 and 0.11 respectively, for the run shown in the Figure. Although the 
temperature profiles, both kinetic and configurational, appear to be well converged, s co id 
increases and Shot decreases, though with significant fluctuations, throughout this relatively 
long run. We turn next to the somewhat simpler appearing version of Nose's idea promoted 
in Reference 36. 



"Weakly Coupled" Logarithmic Thermostat 

About a year ago Campisi et alii [ CZTH ] introduced a logarithmic thermostat much 



like Nose's but without an explicit coupling between the thermostat and the system 
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. So as to apply their idea to the sixty-particle 4 chain we attach a CZTH logarithmic 
thermostat to each of the forty thermostated particles : 

60 40 

Hczth = E(P 2 /2m), + (rf/2M)i + (m/2) In - Si ) 2 + S 2 ) } + $(<?) . 
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FIG. 3: The upper plot shows the configurational temperature profile averaged over times of 3000, 
10000, and 100 000, with the last time shown as open circles. Notice that the configurational tem- 
perature equilibrates well despite the discontinuities in the nearby constrained kinetic temperatures 
shown in the lower plot. The 60 particle kinetic temperatures are represented by small filled circles. 
The open circles in the lower figure are the kinetic temperatures of the Logarithmic-thermostat 
particles. Notice that the "hot" thermostat particles have measured temperatures far below the 
specified T = 1.5. There is no measured heat flux in this simulation. Two billion timesteps, with 
dt = 0.00005, were used to compute these averages. The Hamiltonian is 52.882. 

The total number of differential motion equations to be solved is 200, four for each ther- 
mostated particle ( q, p, s, p s ) and two for each Newtonian particle ( q, p ) . In addition we use 
100 more equations to compute the time integrals of the 60 particle's kinetic temperatures 
and the 40 thermostat variables' kinetic temperatures. The configurational temperatures , 

kT c = ( F 2 )/( V 2 H ) , 

for the 60 particles require 120 additional equations, 60 for ( F 2 ) and another 60 for the 
particles' ( V 2 "H ) . Just as before we set Boltzmann's k and all the particle masses equal 
to unity. We also set the thermostat-variable masses { M } all equal to one. We choose the 
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Time-averaged kinetic and configurational profiles using the Hoover-Leete La- 



grangian/Hamiltonian thermostat to maintain the kinetic energies of the cold and hot regions. 
Despite the successful thermostating of the reservoirs there is no heat flow. The configurational 
temperatures, indicated by open circles, are nearly constant, and are equilibrated with the 20 New- 
tonian particles in the center of this 4> A chain. This simulation includes a billion timesteps with 
dt = 0.001. The value of the Hamiltonian is 49.282. 

"small" parameter 5 = 0.01 in order to avoid the singular behavior of purely-logarithmic 
thermostats. Once again the logarithmic thermostat leads to stiff equations, requiring a 
timestep of 0.00005 for six-figure energy conservation in a two-billion-timestep run. 

With the logarithmic "thermostat" we show separately both the long-time-averaged ki- 
netic and configurational temperatures of the reservoir particles, ( p 2 ) and the kinetic 
temperatures of the corresponding thermostat variables ( p 2 s ). It is disconcerting to learn 
that the temperatures of the thermostat variables and the thermostated variables they are 
assigned to control can differ by a factor of two! Otherwise the basic kinetic-temperature 
results are very like those found in our application of Nose's original Hamiltonian work. 
All of the temperatures are shown in Figure 3. There is no trace of a smooth temperature 
gradient like that generated by the Nose-Hoover motion equations. We conclude that the 
CZTH "thermostat" behaves much like Nose's. The logarithmic thermostats are evidently 
both "stiff" and ineffective. 
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Hoover-Leete Isokinetic Thermostat 



One of us, Bill, met with Tom Leete where Tom worked, at the Aberdeen Proving Ground 
in Maryland. Tom had chosen to work for the Army there because it was "easier" than grad- 
uate school. This was shortly after Leete had finished his 1979 Master's Thesis work at West 
Virginia University, "The Hamiltonian Dynamics of Constrained Lagrangian Systems" . Bill 
was still in the throes of developing the thermostated nonequilibrium shear and heat flow 
treatments he had worked out with his first Ph. D. student William T. Ashurst. Ashurst's 
thesis, "Dense Fluid Shear Viscosity and Thermal Conductivity via Nonequilibrium Molec- 
ular Dynamics", was published in 1974. Bill worked out, with Tom's help, a Hamiltonian 
thermostat which included a nonholonomic constraint on the particle-based kinetic energy, 
K(q)=E(mq 2 /2) : 



^hl( q,p) = 2^jK(p)K(q) - K(q) + ; K(q) constant . 

The resulting Hamiltonian equations of motion allow the "momenta" { p } to change with 
time in the usual way, but constrain the summed up squares of the velocities { q } to remain 
constant, 



£mg 2 = 20kT ; { q = (p/m)yj K(q)/K(p) ; p = F(q) } . 

Generally implementations of nonholonomic constraints lack uniqueness. But a unique im- 
plementation results if one applies Gauss' Principle (of Least Constraint) 17|, |4l|, |42|. Al- 
though the Hoover-Leete motion equations match the Gaussian isokinetic motion equations 
to second order in the time the third-order motions, q, differ. The Hoover-Leete approach 
is a somewhat different way of imposing constant kinetic energy on a Hamiltonian system. 
Dettmann and Morriss derived a Hamiltonian, different to the Hoover-Leete one, which 
satisfies Gauss' Principle and likewise imposes a constraint on the kinetic energy. 

The dynamics from the Hoover-Leete Hamiltonian is straightforward and, applied to our 
heat-flow problem, leads to analogs of the logarithmic-thermostat approaches. As usual, the 
particle masses and Boltzmann's constant are set equal to unity. The motion equations are 
less stiff, with a timestep dt = 0.001 conserving energy with six-digit accuracy for a billion 
timesteps. Just as with the other Hamiltonian approaches to two-temperature mechanics, 
the cold and hot portions of the chain appear to have no influence on their Newtonian 
neighbors. And again, a Hamiltonian approach to thermostating fails, as shown in Figure 4. 
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Configurational Thermostat 



Landau and Lifshitz' expression for configurational temperature 34J, 

kT c = ( F 2 )/( W 2 H > , 

can likewise be applied to the heatflow problem. Figure 6 of Reference 17, for a 600-particle 
<f) A chain, with 200 particles in each of the three regions, shows again that there is no 
tendency for the reservoirs to influence the temperatures of the nearby Newtonian particles. 
We resist carrying out additional simulations for the model here because the equations of 
motion are relatively complicated, requiring, as they do, a Lagrange multiplier to control 
not only T = 0, but also T = . The 600-particle results in Reference 17 lead again to 
the same conclusion: Hamiltonian thermostats are ineffective, unless, as is the case with 
the configurational temperature, the degrees of freedom being thermostated are already at 
the "right" temperature and with the first and second time derivatives of that temperature 
constrained to vanish. 



SUMMARY AND CONCLUSIONS 

Reservoirs based on Nose-Hoover dynamics are an efficient and particularly successful 
route to nonequilibrium properties. The underlying derivation for this dynamics can be 
based on Nose's Hamiltonian, or on the Dettmann-Morriss Hamiltonian, or on the phase- 
space continuity equation necessarily obeyed by any equations of motion. The nonHamilto- 
nian approach provides the simplest route and has been generalized in many useful directions, 



even including Brownian motion 30] . Applied to the 4 model the Nose- Hoover reservoirs 



generate heat flow obeying Fourier's Law, along with realistic temperature profiles. 



It has been argued that heat flow is an unnecessarily demanding thermostat test [43]. 
We strongly disagree. Unless a "thermostat" is capable of transporting heat away from 
"hot" degrees of freedom and transferring heat toward "cold" ones, it is certainly unfit to 
"control" temperature. For this reason we advocate testing any proposed thermostats with 
the highly-useful 4 model investigated in the present work. 

Jones and Leimkuhler 44j have recently studied the usefulness of stochastic forces in deter- 



mining the usefulness (they term this "adaptability" ) of thermostats. They seek thermostats 
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which are both ergodic (in very small systems this is a reasonable request) and which ex- 
actly reproduce their target temperature (we agree, and view this requirement as a necessary 
property of any "thermostat" worthy of the name). Because we believe that determinism 
and time-reversibility are also indispensable thermostat properties, we don't favor stochastic 
thermostat tests. 

Here we have used the 4 model as a test of four Hamiltonian reservoir models, [1] 
Nose's original Hamiltonian, [2] the Logarithmic Hamiltonian reservoir of Campisi et alii, 
[3] the Hoover-Leete isokinetic Hamiltonian reservoir, and [4] the Travis-Braga version of 
Landau and Lifshitz' configurational temperature reservoir. Only the Nose-Hoover reservoirs 
generate a realistic heat flow. The Hamiltonian models are relatively stiff to implement and 
provide no equilibration between the hot and cold reservoirs. Evidently the twin restrictions 
of constant energy (Hamiltonian) and constant temperature so constrain the reservoirs that 
they are unable to influence even their nearest neighbors. The failure of every one of these 
Hamiltonian heat flow simulations underscores the need to generalize mechanics, as did 
Shuichi Nose, in order to treat nonequilibrium problems. Hamiltonian mechanics itself 
cannot provide heat sources or sinks. It is simply too specialized for the realistic treatment 
of nonequilibrium flows. 

It is also noteworthy that the Hamiltonian thermostating methods require one or two 
orders of magnitude more computer time in order to match the accuracy of the Nose- Hoover 
thermostat. 
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